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Introduction 


In Unit 14 on vector algebra we showed how a physical quantity could be 
categorized as either a scalar or a vector, and we spent some time on the algebra 
of vectors, i.e. adding, scaling and multiplying vectors. A scalar quantity has 
magnitude only and can be represented by a real number; for instance, mass and 
temperature are scalars. A vector quantity has two properties associated with it, its 
magnitude and a definite direction in space. We represent a vector on a diagram 
by an arrow or algebraically by components. Velocity and force are examples of 
physical quantities which can be represented by vectors. 


In Unit 14 the values of the scalars and the components of the vectors did not 
depend on where they were measured. 


In this unit we shall consider scalars and vectors which have different values at 
different points in a region of space. For instance, you know from practical 
experience that the temperature in a bowl of hot soup varies quite significantly 
from the middle to the edge of the bowl. That is why we try to drink the soup 
from the edge rather than from the middle where we are likely to get a burnt 
mouth! Temperature is a scalar quantity and its value depends on where in the 
soup it is measured. The temperature distribution is an example of a scalar field. 


Similarly, the force on a nail in the region surrounding a magnet will depend on 
how near to the magnet we bring the nail. The magnetic force field of the magnet 
has a definite magnitude and direction which varies with the position of the point 
at which it is measured. This is an example of a vector field. 


This unit is about scalar and vector fields. 


Study guide 


This unit is divided into four sections and they should be studied in the order in 
which they appear. The audio-tape is part of Subsection 3.3 and the television 
programme should be viewed at the beginning of Section 4. If you have not had 
time to study the first three sections before watching the television programme 
then a quick skim through Sections 1 and 3 should be sufficient for you to follow 
and enjoy the programme; then return and study the unit in more detail. 


In Section 1 we define scalar and vector fields and show how we can describe 
them geometrically. In Section 2 we define the spatial rate of change of a scalar 
field called the gradient of that scalar field. Section 3 introduces a new kind of 
integral of a vector field called the line integral. You should listen to the tape for 
this unit while studying this section; it takes you through several line integrals to 
help you gain confidence at evaluating them. We then, in Section 4, define a new 
vector field related to the spatial rate of change of the old vector field and called 
its curl. This definition uses the line integral introduced in Section 3. The television 
programme, which is associated with this section, shows you the physical 
significance of curl by considering the velocity vector fields of three fluid flows. 


Finally, a note to those students who have previously studied MS283. You should 
be familiar with the material contained in this unit. You are recommended to 


work through the exercises and to revise those sections whose exercises you find 
difficult. 
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1 Scalar and vector fields 


1.1 A scalar field defined 


Figure 1 shows the temperature at different points of a flat plate when unevenly 
heated and cooled. The value of the temperature depends on where on the plate it 
is measured. If we define Cartesian coordinate axes Ox and Oy as shown, then the 
variable T representing the temperature will be a function of x and y; ie. 

T= T(x, y). This function of two variables is just an example of the concept of a 
function of more than one variable introduced in Unit 25. Figure 1 


KR 
heat source 


As another example, Figure 2 shows the (atmospheric) pressure distribution over 
the U.K. at midday on 30 September 1980. 


‘004 


Figure 2 


The solid lines are called isobars and join up places where the atmospheric 

pressure is the same. Atmospheric pressure at a given time is a scalar quantity and 

its value depends on where in the U.K. it is measured. For instance, at midday on 

30 September 1980, in Manchester it was 1020 millibars and in Penzance it was 

1028 millibars. 1 millibar = 100 Nm? 


If we choose London as an origin and set up a Cartesian coordinate system with 
Ox pointing East and Oy pointing North, then each place in the U.K. can be 
represented by its coordinates (x,y). Then the variable P representing the 
atmospheric pressure is a function of two variables x and y; i.e. P = P(x, y). 


These two examples are pictorial representations of scalar fields. In each of the 
examples above we have: (i) a scalar quantity (e.g. temperature or atmospheric 
pressure), (ii) a region of space over which each scalar is defined (e.g. the flat plate 
or the U.K.), and (iii) the values of the scalar quantity at different points of the 
region. 


In each of the above cases, the variable representing the scalar quantity is a 
function of the two variables x and y; these are two-dimensional scalar fields. The 
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temperature distribution in the room where you are reading this unit would be an 
example of a three-dimensional scalar field, which we now define: 


If @ is a function which associates a unique scalar with each point 
in a given region of space, then @ is called a scalar field function 
or simply a scalar field. 


The domain of the function @ is the given region of space and the codomain is the 
set of all scalars. 


Example 1 


The light intensity from a point source of light, such as a small light bulb, is 
inversely proportional to the square of its distance from the bulb. The intensity at 
distance 1 metre is I). Specify a scalar field function which models the light 
intensity and state the region of space over which it is defined. 


Solution | 
Suppose that we represent the light source by a point at the origin, and the light 


s e ° ° ° 1 
intensity at a point distance r metres from the light source by J; then J oc — and 
: r 


I 
I = Ip when r = 1. Thus I = 5. 


Io 


The scalar field function is J = 3 (o in Cartesian coordinates). The 


0 
: : (x? + y? + 27) 
function is not defined when r = 0 and so the region on which the scalar field is 
defined is the whole of space with the origin (0,0,0) removed. Note that this is a 


three-dimensional scalar field. 


Exercise 1 


A thin circular sheet of metal of radius 4 metres has a concentric hole of radius 1 metre cut 
out of it (see Figure 3). The inner perimeter is maintained at a temperature of 0°C and the 
outer perimeter is maintained at a temperature of 100°C. At a point P on the sheet, given 
by position vector r, the temperature is proportional to log, r, where r = |r|. Specify the 
scalar field function which models the temperature and state the region of the plane over 
which it is defined. 


[Solution on p. 44] 


1.2 Contour curves and contour surfaces 


The scalar fields in the last subsection were presented in three different ways. For 
the first scalar field of temperature over a plate, as given in Figure 1, we 
represented the data as values of the temperature at a few points on the plate. For 
the second scalar field of atmospheric pressure in the U.K., as given in Figure 2, 
we represented the data by isobars (or lines joining points at which the 
atmospheric pressure is the same) on a map of the U.K. (In both cases this is not 
a complete representation of the scalar field because for instance, there is a 

_ temperature associated with every point of the flat plate, and not just at the points 
shown in Figure 1.) Finally, in Example 1 and Exercise 1 the scalar fields were 
given by a formula of some type. 


Of these modes of representation, the isobar example gives a very good way of 
visualizing pictorially the physical situation. Figure 4 shows the temperature 
distribution over an unevenly heated flat plate similar to that shown in Figure 1. 
However instead of giving the value of the temperature at different points as we 
did before, we have drawn curves through those points on the plate where the 
temperature is the same. Of course we cannot show all possible curves because 
that would require a curve through every point on the plate, so we choose certain 
values of the temperature and the curves give us a pictorial idea of the scalar field. 
So, for example, for each point on the curve AB the temperature equals 300°C. 


We can define a two- 
dimensional scalar field 
similarly, but then the domain 
is a region of the plane and 
not a region of space. 


Figure 3 


heat source 
Figure 4 
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For a temperature field the curves are called isotherms, and in this example we 
have drawn the isotherms for multiples of 50°C. Notice that in one part of the 
plate the isotherms are much closer together than for the rest of the plate; the 
temperature is changing more rapidly in the region where the isotherms are close 
together. 


We now extend this pictorial representation for pressure and temperature to more 

general two-dimensional scalar fields. Suppose we have a two-dimensional scalar _ 

field (x, y) defined over a region R. We can represent this scalar field pictorially \ 
by drawing curves through those points of R which have equal values of @. Figure 

5, for example, shows four such curves where C,, C,, C3 and C, are constants. 

The curves that we have drawn are called contour curves. Figure 5 


For example, consider the scalar field 
@=xX +¥ 


defined over the rectangle —1 < x < 1, —2 < y <2. The contour curves are 
defined by 


@ = constant, or x? + y* = constant. 


These curves are circles or parts of circles drawn within the rectangular region for 
which @¢ is defined. Figure 6 shows the contour curves 


&(x, y) = 0.25; O(x%y)=1;  o(% y)= 2.25; (x, y) = 4. 


Again we cannot draw all possible contour curves, only some of them. 


Exercise 2 


Draw some of the contour curves to represent pictorially the scalar field ¢ = xy defined 


over the disc x? + y? < 4. 
Figure 6 
Solution on p. 44] 


You are probably familiar with the word ‘contour’ from Ordnance Survey maps. 
These are lines on a map joining up places which are at the same height above sea 
level. The contours can give us a lot of information about the lie of the land. Here 
is an example; see if you can answer Exercise 3, just from the contour curves. 


Figure 7 


Exercise 3 
(i) Where would you expect to find (a) a hill, and (b) a lake? 


(ii) Where do you think that the land is fairly level? 


(ii) Suppose you set off from point A and walked in the direction of the arrow, would you 
expect an easy walk or a hard climb? 


(Solution on p. 44] 


This exercise illustrates one important feature of a scalar field which can be 
obtained from its contour curves. The rate of increase (or decrease) of a scalar 
field can be deduced from the closeness of the contour curves. If the contour 
curves are close together then the scalar field is changing quite rapidly in that part 
of the region. 
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So far in this subsection we have considered two-dimensional fields. Now we turn 
our attention to three-dimensional scalar fields. A pictorial representation of a 
two-dimensional scalar field (x, y) is given by the contour curves defined by 

(x, y) = constant. We can represent a three-dimensional scalar field (x, y,z) in an 
analogous way. If we join up points for which the three-dimensional scalar field 
(x, y,z) has the same value, we obtain the surface with equation $(x, y,z) = C 
where C is constant. The family of surfaces given by $(x, y,z) = C, for different 
values of the constant C, are called the contour surfaces of the scalar field and they 


give a pictorial representation of the scalar field. 


For example, the contour surfaces of the scalar field I = I)/(x? + y? + z’) given in 
Example 1 are concentric spheres with centres at the origin. Two of the contour 


x 


surfaces, for J = 4 and J = 9, are shown in Figure 8. Figure 8 

1.3. A vector field defined <a = 
x a oa ae 

Figure 9 shows what happens to iron filings when placed near a bar magnet. The --===sjS___—sNE==---- 

iron filings have aligned themselves in a symmetrical pattern and they are . S = ee _ 

experiencing a force due to magnetism. If we place small compasses at various ee Se 

points near the magnet we find the arrow on a compass orientates itself along the Figure 9 

iron filings. The compasses show us the direction of the magnetic field at each 

point. The magnitude of the field decreases as we go further away from the See 

magnet. The arrows in Figure 10 represent the magnetic field due to the bar n 

magnet at various points near to the magnet. At each point in a region close to Saree 


the magnet there is a unique vector called the magnetic field at that point. This is 
an example of a vector field. 


Here is another example. The arrows in Figure 11 represent the surface velocity at ? 
various points in a river. aa as 
—s 


Figure 10 


Figure 11 


Near the river bank the water is almost at rest, so the arrows are somewhat 
shorter in length (indicating a smaller velocity magnitude near the edge of the 
river). The velocity at a particular point is the velocity that a small floating object, 
such as a leaf, would have at that point. At each point on the river’s surface there 
is a unique velocity; this is an example of a velocity vector field. 


You have seen something like this before in Unit 2 where you met direction fields. 
These were a distribution of orientations in space, i.e. directions without a sense 
being attached, but no magnitudes were associated with the elements of the 
direction field. We now have a distribution of vectors in a region of space, 1.e. 
quantities with magnitude and direction, which may take on different values at each 
point in the region. This leads to the following definition: 


If F is a function which associates a unique vector with every point 
in a given region of space, then F is called a vector field function 
or simply a vector field. 
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The domain of the function is the given region of space and the codomain is the 
set of all vectors. To describe a vector field completely we require the following: 

(i) the region of space on which the vector field is defined, and (11) the magnitude 
and direction of the vector at each point of the region. 


Example 2 

The Earth’s gravitational field can be considered as an example of a vector field: 
as you will see in Unit 30, the force on a body which is at a distance r from the 
centre of the Earth has magnitude inversely proportional to the square of the 
distance r and is directed towards the Earth’s centre. Specify the vector field and 
the region of space over which it is defined. 


Solution 


The vector quantity is the gravitational force acting on the body due to the Earth. 
This force acts at points (x, y,z) in space such that 


J/x*+y*+27>R 


where R is the radius of the Earth, and x, y, z are Cartesian coordinates with 
origin at the centre of the Earth. The vector field function is — F,e, where e, is a 
unit vector whose direction is from O to P (see Figure 12). The magnitude F, is 
proportional to 1/(OP)?, so 


= Ss 
* OF —i oy +2 } Figure 12 


x 


for some constant C, and the direction of F is along the vector 
YP ° . 
PO = —(xi+ yj + zk) 

so that 


eee 


—_ Pea —C (xi + yj + zk) 
‘loP| (x? +y +27) fry ge 
= Coe + yt Zh) 
(x? ae y? e y ‘ 
The region over which the vector field F is defined is the set of points (x, y, z) 
outside the sphere x? + y? + z? = R?. 


This is an example of a three-dimensional vector field because the function F is 
defined on a three-dimensional region of space. 


The components of the vector field function relative to a given Cartesian 
coordinate system are themselves scalar field functions defined over the same 
region of space as F. Thus if F is a three-dimensional vector field function defined 
on a region D and 

F(x, y, Z) = F; (x, y, z)i + F,(x, y, Z)j a2 Pix, z)k 
then F; (x,y,z), F:(x, y,z) and F;(x, y,z) are scalar field functions also defined on D. 


The following exercise gives an example of a two-dimensional vector field. 


Exercise 4 


A flat horizontal metal circular plate of radius 10cm has a small hole in the middle through 
which is injected water at a constant rate. It is found that the velocity of the water at any 
point P on the plate is directed outwards along the line joining P to the centre of the plate 
O and has magnitude inversely proportional to the distance from O. Specify the vector field 
and the region over which it is defined, in terms of Cartesian coordinates with origin at O. 


[Solution on p. 44 | 


1.4 Vector field lines 


The arrows on the diagrams in the last subsection are one way of describing a 
vector field pictorially. The pictures are however incomplete because there is a 
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vector quantity associated with every point in the region. A more complete picture 
can be obtained by drawing continuous curves in the region such that at any point 
the tangent to the curve is parallel to the direction of the vector field function at 
that point. These curves are called the (vector) field lines of a vector field. For 
instance, the orientation of iron filings near a bar magnet suggests a family of 
continuous curves. Two such curves are shown in Figure 13. At any point on one 
_of the curves the orientation of a compass needle would be along the tangent to 
‘the curve; one is shown in the figure. 


Exercise 5 
Sketch the field lines for the Earth’s gravitational force field. 


[Solution on p. 44 | 


When the vector field represents a force the field lines are often called lines of 
force. 


Here is another example of field lines. The arrows in Figure 14 represent the 
velocity of the surface of.a river at various points. If we drop some leaves onto the 
river they will float. If we draw a line showing the path of one of these leaves, then 
at each point on the path, the tangent to the path is parallel to the direction of the 
velocity vector. This line is a field line for the velocity vector field. For the motion 
of a liquid the field lines are often called flow lines. 


You have met field lines before in Unit 2 of this course. The trajectories of a 
direction field are analogous to the field lines of a vector field. 


The field: lines give a convenient way of describing a vector field pictorially, and a 
diagram showing the field lines often helps us to describe the vector field which 
they illustrate. However the field lines tell us only the direction of the vector field 
and not the magnitude. 


Example 3 
Describe briefly the vector field corresponding to the field lines in Figure 15. 


Solution 

The vector at any point P is along the tangent to the circle through that point (see 
Figure 16), and this direction is perpendicular to the vector OP. Since an 
anticlockwise sense around the circles is indicated by the arrowheads, the direction 
of F is that of —sinOi + cos6j. If P is the point (x, y) then sin® = y/./x? + y* and 


cosO = x/,/x? + y*, so the direction of F is that of (— yi + xj). The vector field 
function has the form F = f(x, y)(—yi+ xj) for some function f(x, y) taking 
_ positive values only. 


Exercise 6 
Describe briefly the vector field corresponding to the field lines in Figure 17. 


[Solution on p. 44] 


Summary of Section 1 


1. If @ is a function which associates a unique scalar to each point in a given 
region of space, then ¢ is called a scalar field function or a scalar field. 


2. A scalar field in two dimensions is defined over a region of the plane and can 
be represented by a function of two variables, #(x, y). Such a scalar field can 
be described pictorially by contour curves; these are curves drawn on the x, y- 
plane with equations given by (x, y) = constant. 


3. Ascalar field in three dimensions can be represented by a function of three 
variables, (x, y,z), and can be described pictorially by contour surfaces. These 
are surfaces with equations given by (x, y, z) = constant. 


4. IfF is a function which associates a unique vector with every point in a given 
region of space, then F is called a vector field function or simply a vector field. 


aE 


Figure 13 


Figure 16 


aK 


Figure 17 
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5. A vector field can be described pictorially by field lines. A field line is a curve 
whose tangent at each point is in the direction of the vector at the point. 


2 Differentiating scalar fields 


In the study of functions of more than one variable in the last unit, certain 
properties of the function were best understood in terms of the partial derivatives 
of the function. For instance, at a stationary point of a surface the first-order 
partial derivatives are zero. Whether the stationary point is a local maximum, 
local minimum or saddle point depends on some combination of the second-order 
partial derivatives. 


Certain characteristics of scalar and vector fields can be described in terms of 

the spatial rate of change of the fields, and we shall see that these involve a 
combination of the first-order partial derivatives of the field functions. The 
properties that we are going to introduce in this and the next section provide us 
with relationships between a scalar field and a vector field, and between different 
vector fields. These are fundamental properties of fields that will be used in the last 
section and in Unit 30 to characterize and describe the influences of various 
physical quantities in the real world. 


2.1 The gradient of a scalar field 


A scalar field is described pictorially by contours, i.e. contour curves for a two- 
dimensional field and contour surfaces for a three-dimensional field. The 
closeness of the contours tells us quite a lot about how the field is changing in a 
particular direction. For instance, land elevations are represented by contour 
curves drawn on a map. If the contour curves are close together then there is a 
steep hill in that region and if the contours are well separated then the land is 
reasonably flat. 


The change in a scalar field in different directions can also be described by looking 
at the contours. Suppose we consider the elevation of the land above sea level, 
represented by the function @(x, y), as an example of a scalar field. Figure 1 shows 
a hill and its representation on a map by contour curves. 


p 


a 
a 


x 


Figure 1 


If we travel round a contour curve from A to B say, then the value of the function 
(x, y) does not change because a contour curve is defined by (x, y) = constant. 
Now suppose that we put a ball at the point A on the hill. Experience tells us that 
the ball will begin to roll down the hill along the line of steepest slope, in other 
words the direction in which the scalar field is changing most rapidly. We would 
probably expect that this direction is perpendicular to the tangent to the curve AB 
at the point A. The acceleration of the ball will depend on how steep the slope of 
the hill is near the ball’s position. At each point of the hill the magnitude and 
direction of the initial acceleration of a ball defines a vector field. We would expect 
this vector field to depend on how the scalar field of land elevations changes from 
point to point. It is this change in a scalar field that is described by the gradient. 


By ‘spatial’ we mean with 
respect to space coordinates; 
that is, not with respect to 
time. 
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In this subsection we are going to make this idea more precise. We begin by 
adopting the following definition of the gradient: 


At each point in the domain of a two-dimensional scalar field ¢ 
we can evaluate the vector 


giving a vector field called the gradient of the scalar field @ and 
denoted by grad @. 


We will see that the gradient of a scalar field tells us how the scalar field is 
changing. But before we do that, the following examples and exercises show that 
finding the value of grad @ at a point is relatively straightforward. 


Example 1 
Evaluate grad @ at the point (0,1) when @ = log,(x + 2y). 


Solution 
The vector field grad ¢ is formed by evaluating the partial derivatives of @ of first 
order: 
d = log.(x + 2y); 

0p 1 0 2 

2 a and Bh = . 

Ox x+2y Ody x+2y 

ee? ‘aoe © 

So grad @ = F z mL +- F : hi 


Putting x = 0 and y = 1 we have, at the point (0,1), grad @ = 41 + j. 


Exercise 1 
Evaluate grad ¢ at the point (— 1,2) when ¢ = x’y. 


[Solution on p. 44 | 


Example 2 
Evaluate grad ¢ at the point (1,2) when @ = x? + y*. Show that the direction of | 
grad ¢ at (1,2) is normal to the contour curve through the point (1, 2). oo means ‘at right angles 
Solution 
- The vector grad ¢ at (1,2) is found by evaluating the partial derivatives of @ of 
first order at the point (1,2): 
GSI, 

og 0p 

a and an at (1,2); 

op og 

oy = 2y and éy = 4at it 2). y 
So at (1,2) P(1,2) 

0p. , Op 
d@ = —i+ —j = 21+ 4j. 
grad @ a ae 1+ 4) 


Let P be the point (1,2). The contour curve through P (see Figure 2) is a circle of 
radius ./5. The normal to the circle through P is in the direction of the vector 


OP =i+ 2j = $(2i + 4j), so at (1,2) grad @ is in the same direction as this normal 
vector. Figure 2 
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In fact we will now show that the direction of the gradient of any scalar field ¢ at 
a point is along the normal to the contour through that point. 


Suppose that ¢ (x, y) is a scalar field in two dimensions and consider one contour 
defined by ¢ (x,y) = C where C is some constant, and suppose that P is a point 
on the curve with coordinates (xo, Vo) (see Figure 3). 


In Unit 14 we introduced the idea of a curve in the plane as the representation of 
some vector function of a scalar t; so the position of a point on the contour ¢ = C 
can be written as 

r=x(t)i+ y(t)j 


where t is a scalar, and as t varies the point moves along the contour. A tangent 
vector to the contour curve is then 


dr _ ax. dy 
dt "die det 


With this representation the value of the scalar field (x, y) as we travel along the 
contour curve is also a function of t. Using the chain rule of Unit 25 we can write 


dp _ op dx dgdy 
dt 0x dt Oy dt 


which is the dot product of two vectors, i.e. 


db _ a : “;). (i 4 2 dr 


dt Ox dy! de} eat 
where ae ae 
= a = grad @¢ 


d d 
Since @ has constant value along the contour, © = 0; hence g = = (0, and the 


dr pe 
two vectors g and - are perpendicular. Thus the direction of grad @ is 
perpendicular to the tangent vector, i.e. grad @ is along the normal to the contour. 


Exercise 2 
Find a vector whose direction is normal to the curve x? — 2xy + y? = 9 at the point (0,3). 


[Solution on p. 44] 


For a three-dimensional scalar field represented by the field function (x, y, z), the 
gradient of @ is defined by 


Example 3 C 
Consider the scalar field (x, y, Zz) = 02+ y+ 242 
where C is constant. Show that the gravitational field F = — F,e, defined in 


Example 2 of Section 1 can be expressed in terms of this scalar field ¢ by the 
relation F = grad @¢. 


Solution 
C 0d Cx 
ee ee 
0p Cy Ty cs 


moo and —> = — 


oy (x? + y? + 27)3/2 az (x? + y? + 22)3/2 ’ 


dr 
dt 


P(X ,¥o) 


contour 


o=C 


Figure 3 
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so that 


C <i 
grad d = Zs? (xi+ yj + zk). 


168 + 42 
This is the vector field function describing the gravitational field at the point 
(x, y, Z). 


This example illustrates an important application of the gradient. Scalar and 
vector field functions that we use to model physical quantities in the real world are 
often related in this way; i.e. some vector field function is just a scalar multiple of 
the gradient of some scalar field function. 


For instance, the gradient of the temperature field on the unevenly heated plate in 
Section 1 is a vector field which determines the flow of heat across the plate. In 
Figure 4 the dotted lines represent the contour curves and the arrows represent 
the gradient vector at various points on the plate. 


Exercise 3 
Evaluate grad @¢ at the point (—1,1,0) for @ = x*y + 27x. 
[Solution on p. 44 | 


At a point P on a contour surface describing a three-dimensional scalar field, the 
gradient is perpendicular to the tangent plane at P. 


Example 4 — 


Find the direction of a normal vector to the surface n;x + n,y + n3z = C at the 
point (1,0, 2). 


Solution 


The surface n,x + n,y + n3z = C is a contour surface of the scalar field 
p(x, y,z) = nx + ny + n3z and the gradient of ¢@ evaluated at the point (1,0, 2) is 
normal to the contour surface. We have 


ad ao 


—=n,, — =n, and 
Ox 


Op _ 
dz 
so that grad @ = n,i + n.j + n3k. This vector is the same for every point on the 


surface. This result should not surprise you because the equation 
nyX + ny + n3z = C is the equation of a plane, as we saw in Unit 14. 


N3 at (1,0, 2}, 


The gradient of a scalar field ¢ is defined in terms of the partial derivatives of ¢ of 
first order and its direction is along the normal to the contour curves for a two- 
dimensional scalar field or along the normal to the tangent plane for a three- 
dimensional scalar field. So in some sense it is telling us about the change in ¢ as 
we go along the direction of the normal. 


Consider again a two-dimensional scalar field 6; we now investigate the rate of 
change of @ in the direction of an arbitrary unit vector e. 


Figure 5 shows several contour curves for this scalar field and at a point P a unit 
vector e. Suppose that the curve C is the path of a particle moving through the 
domain of the scalar field with unit speed, and that the direction of the tangent to 
the curve at the point P is given by the unit vector e. Then the velocity of the 
particle as it passes through P is equal to e. Using the idea of a curve as the 
representation of a vector function of time t as introduced in Unit 14, the position 
of the moving particle on the curve can be written as 


r= x(t)i+ y(t)j 


and then e = si 
dt 


heated corner 
Figure 4 


@ = constant 


(contours) 
ee 


Figure 5 
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Using the chain rule of Unit 25 we have the rate of change of ¢ with respect to t, 
at the point occupied by the moving particle, as 


dd 0p dx 0p dy 
dt dxdt dydt = eter | 


d 
Now grad @ is along the normal to the contour at P, and sak 


- is the velocity of the 


moving particle which as we have seen is equal to e. 


Suppose that @ is the angle between grad @ and e (see Figure 6). Then the above 
equation becomes 


7 aes (grad p)-e = |grad ¢||e| cos 6 = |grad ¢| cos @ 


since e is a unit vector. 


Considering how this expression varies with @ we see that the rate of change of ¢ 
is zero when 0 = 7/2, i.e. along the contour. This is what we would expect because 
the contour is defined by ¢ = constant. The maximum rate of change in @ occurs 
when 6 = 0, i.e. along the direction of the vector grad ¢, and this maximum rate 
of change of ¢ is |grad ¢|. This means that the gradient of a scalar field is the 
vector field which represents the magnitude and direction of the maximum (spatial) 
rate of change of the scalar field. This maximum rate of change is along the normal 
to (i) the contour curve for a two-dimensional scalar field and (ii) the contour 
surface for a three-dimensional scalar field. 


Exercise 4 
For @ = x’y?z? find 


(i) the maximum rate of change of ¢ at the point (—1,1,1), and the direction in which 
this maximum occurs; 


(ii) the rate of change of ¢ at the point (2,1, —1) in the direction of the unit vector 
e=2i+ $k. 
[Solution on p. 45 | 


2.2 The vector operator V 


We now introduce a special notation for the gradient of a scalar field. Instead of 
grad @ we often write V@ where the symbol V is pronounced ‘del’ or ‘nabla’ and 


é 


0 0 
V stands for i— +j—+k—. 
oy Oz 


Ox 


When written to the left of the symbol representing the scalar field function 
(x, y,z) it means ‘find the first-order partial derivatives of @ and form the vector 


¥6= i+) +a 


The outcome of this operation is the gradient of the scalar field function. Usually 
we call V@ not ‘del d’ but ‘grad @’. V is not a vector by itself but is called a 
differential vector operator implying that it needs something to ‘operate on and 
differentiate’ in a special way. 


Exercise 5 
We know that if f and g are two functions of a real variable t then 


“= Sat fe 


What would be the analogous result for V(¢w) where ¢ and w are two scalar fields? 
[Solution on p. 45 | 


grad d 


Figure 6 


@ = constant 
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Summary of Section 2 


eo eee, | 
1. The vector field ae + ae jt ok formed from a scalar field ¢ is called the 
x y Z 


gradient of @ and this is written as grad @ or V@¢. 
2. Properties of grad ¢: 


(i) the direction of grad ¢@ at a point P in the domain of the scalar field @, is 
either along the normal to the contour curve drawn through P for a two- 
dimensional scalar field or along the normal to the tangent plane of the 
contour surface drawn through P for a three-dimensional scalar field 


(11) the magnitude of grad @ at a point P is the maximum rate of change of @ 
at the point P. 


3. The rate of change of a scalar field ¢ at a point P in the direction of a unit 
vector e is given by (grad ¢)-e. 


3. The scalar line integral 


3.1 What is the scalar line integral? 
You are familiar with the idea of integrating a continuous function over an 


b 
interval [a,b] of the real line, i.e. you can evaluate integrals of the form f (x)dx 


which are called definite integrals. The definite integral is defined as the limit of a 
sum in the following way: 


We divide the interval [a,b] into N subintervals and in each subinterval we select 
a point: 
x, in the first subinterval, 


x, in the second subinterval, 


x; in the ith subinterval 


and so on. 
Then we form the sum 

F(x )Ox, + i g)ehn + + FO + + fen Xn 
where 6x; is the width of the ith subinterval. 


The limit of this sum as we increase the value of N indefinitely is defined to be the 
definite integral of f from a to b and is written 


[ reoay = lim y I (x;)0x;. 


This section is about generalizing the definite integral to give a rather special 
integral of a vector field function that occurs in applied mathematics; it is called 
the scalar line integral. You will see one application of it in Section 4 of this unit. 


Consider a vector field F and let A and B be two points in the domain of F joined 
by a curve C (see Figure 1). Divide the curve into n segments and suppose P, and 
P.,, are two points on the curve with position vectors r; and r; + 6r;, where PP, , 
is one segment. At each point on the curve C the vector field has a particular 
magnitude and direction; at P, let the vector field function be F;. The scalar line 


integral of F along C is defined to be lim )° (F;-6r;), and denoted by 


i=1 


A Fay, or sometimes by | F- dr. 
C 


AB 


The definite integral is defined 
in M101, Block III; MS283, 
Block III; TM281, Block IV. See 
also the next unit in this 

course. 


Figure 1 


If F = f(x)i and C is the 
interval [a,b] of the x-axis 
then 


[ Fear —— f food 
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The curve C is called the path of the integral. The path of a scalar line integral is 
said to be an open path if the two end points do not coincide, as in Figure 2, or a 
closed path if the two end points do coincide, as in Figure 3. 


We can integrate along the closed path in two directions, ice. Starting at A and C 
ending at B or starting at B and ending at A. We usually adopt the convention of 
going anticlockwise round closed curves and to remind us of this we write the 


scalar line integral along a closed curve C as $ F dr. A 
Figure 2: an open path 

The scalar line integral has many applications. For instance, if a vector field and a 

scalar field @ are such that F = V@, then given F we can find @ by evaluating the 


scalar line integral | F: dr; we shall be more precise about this later in this 
C 


section. In Section 4 we shall describe one of the properties of a vector field, called 


the curl, in terms of a scalar line integral. 3 : 
Foe Figure 3: a closed path 
Before we see how to evaluate scalar line integrals we introduce a method of 


representing every point on a path in terms of a single scalar parameter. 


3.2 Parametrization of curves 


In two dimensions a curve can often be described by an implicit equation 

f(x,y) = 0. For example, f(x, y) = y — 2x — 3 = 0 is the equation of a straight line 
or f(x,y) = x? + y? — 9 =0 is the equation of a circle of radius 3 and centre the 
origin. In three dimensions this method is not suitable. An implicit equation 

f(x, y,z) = 0 describes a surface, not a curve. 


We shall adopt a method to represent the paths of line integrals which applies to 
curves in two and three dimensions; the method is based on an idea introduced in 
Unit 14. 


In Unit 14 a vector function f of one real variable (t, say) was described pictorially 
by a curve in space so that any point on the curve had a position vector 


r= xi+ yj + zk = f(1), 
or in terms of the components, 
f(t) = fA(Oi+ AOj+ AOk, 
giving the coordinate equations 
x = f;(f), y = f,(t), z= f3(t). 


These coordinate equations are called the parametric equations of the curve, and t 
is called the parameter ; put in a value of t and we have the position vector of a 
point on the curve. 


A parameter familiar to you from mechanics is time. As a particle moves in space 
subject to certain forces, its position vector is a function of time and the Cartesian 
coordinates of the particle’s position at any time serve as the parametric equations 
of the path along which the particle moves. For instance, in a problem of a 
projectile moving in a vertical plane, the position vector of the particle at time t 
after projection is 

r= (uot i Xo)I “> (—4et? = a Vol + Zo)k 
where (Xoi + Zok) and (uoi + vgk) are the initial position and initial velocity of the 
particle respectively. The path described by the particle has parametric equations 

X =Uol + Xo,. y = 0. z= 4Fgt* + vot + 2 


and these equations represent a parabola in the x, z-plane. However, the parameter 
that we use in general need not have any physical significance. 


Example 1 
Consider the parametric equations 


x=2t, y=4t?+1, z=2  (t any real number). 


18 | MST204 26.3 


Find the coordinates of the points for t= —1, —4, 0,4, 1, and sketch the curve 
that these parametric equations represent. 


Solution 


The coordinates of the points are found by substituting each value of t in turn into 
the equations giving the points (—2,5,2), (—1,2,2), (0,1,2), (1,2, 2), (2,5,2). Every 
point on the curve has z = 2 so that the curve lies completely in the plane z = 2, 
as shown in Figure 4. 


Figure 4 


The parametric equations represent the parabola y = x? + 1 drawn in the plane 
z= 2. 


Exercise 1 
Consider the parametric equations 


x=acost, y=asint, z= 0, (0<t< 2n). 


Find the coordinates of the points for which t = 0, t = 2/4, 2/2, 32/4, x, 3n/2, and sketch 
the curve that these parametric equations represent. 


[Solution on p. 45 | ? 


The parametrization of a curve is not unique. For instance we can show that the 


points on a circle of radius a and centre the origin in the x,y-plane can be defined 
by: ie 


(i) x = acost, y=asint (0 <t < 2z) 


2at _ afl —t) 


(ii) x = ——, Ps 
- sa 1 +t? y 1+ t? 


(t any real number). 

In the above examples the parametrization of the curve has been given; but in 
two dimensions a path may be given as an implicit equation and finding a 
parametrization is not always easy. One method that will sometimes work is to 
let x = t (or y = t) and then solve the implicit equation f(t, y) = 0 (or f(x, t) = 0) 
for y (or x) as a function of t. 


For example, the general equation of'a straight line is y — ax — b = 0 and if we let 
x = t then y = at + b. These are the parametric equations of the straight line. 


However, this method does not always give a parametrization of the. entire curve. 
For instance, for the circle x? + y* — a” = 0, if we let x = t, (for —a <t < a) then 


y* =a’ — ¢’. Solving for y gives y = +./a? — t*. The solution y = +./a? — t? 
gives a parametrization for the upper semicircle (y > 0) and y = — ./a? — t? gives 


a parametrization for the lower semicircle (y < 0), but the method does not give a 
parametrization for the whole circle. The parametrization using the trigonometric 
functions, as in Exercise 1, does parametrize the whole circle. 
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Example 2 
The general equation of an ellipse with centre the origin in the x,y-plane is 
. 2 
“z+ oy 1 = 0. Find a parametrization for this curve. 
Solution 


We could try letting x = at (for —1 <t <1) so that y* = b?(1 — t?), but then we 
would have to treat the two halves of the ellipse for y > 0 and y < 0 separately. 
For y > 0 we have y = b,/1 — t? and for y > 0 we have y = —b,/1 — #?. 


A parametrization which covers the whole curve is to let x = acost (for 
0 <t <2n), and then y = bsint. (There are other parametrizations as well.) 


What we are going to do is to use the parametrization of a curve to help us in 
the evaluation of a scalar line integral. Any parametrization will do, but for some 
the resulting integration is much easier than for others. 


3.3 Evaluating the scalar line integral (Tape Subsection) 
On the audio-tape that follows, we shall'see how to evaluate the scalar line integral, B 


| F - dr, for different vector fields F and different paths C. But first we shall show 
Cc 


that by representing the path in terms of a parameter t, the scalar line integral can 
be written as a familiar looking integral in terms of t. 


We defined the scalar line integral as the limit of a sum of the terms F;- or; as the 
number of line segments gets larger and larger. The position vector r; of the point 
P. on the path C (see Figure 5) can be given in terms of a parameter t by Figure 5 


r; = r(t;) = x(t;)i + y(ti)j + 2(t,)k, 
where t; is the value of t corresponding to the point P. 


As we move along the curve, t changes in value and the change in r; for a small 
change in t¢ is just 


dx; dy; dz; dx; dx 
or; = — d6ti + —dtj + — otk By —- we mean — evaluated 
dt i es ic 
at t;, that is, evaluated at the 
. d 
and we can write point P;, and similarly for = 
> F;°6 ae ee ame ot. and dz, 
i=1 i=1 dt 


In the limit as n increases indefinitely we can write 


wo (eae) = [a4 


where t, and t, are the values of t at the points A and B respectively. 
The integral on the right hand side can be converted to a familiar looking integral: 


d 
if we let Fs = f(t) then the scalar line integral becomes 


| f(t) dt. 


The audio-tape develops an algorithm for evaluating scalar line integrals; it shows 
how to make the transformation from 


t2 
[ Fear to | [r-S} ae, 
C : ty dt 


Start the tape when you are ready 
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2 A Tiso- dimensional fine wntegral 


Evaluate { Por 
| e 


where F = (2x+ yt | 


and ¢, 1s the straight line belween (2,0) and (0,2) 


C2 Introduce € 


= ax +b 


b = y intercept = oe 


(2,0) * 


equation of line 13 y= 
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gradient oo =| 
y infercep! b= 2 


equation y= 2-x 


4) Parametrize 


y=2-x 

x=C€ 

y= 

=e goes From (2,0) To (0,2) 


C goes From to ee 
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o Substitution 


im 
H 


my 

it 

% 
Se 
oo 

<< 


Faay (= 2 ond & = © 


i 7 d 
wee Mewar = [da 
a 
Solution al the end of the Frames 
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: © Evaluation 


[@ + 2t) dt 


2 


a Line integrals 


look aT path 


Step | : parametrize of integration 


write all the 


substitute variables x,y (etc) 
in terms of parameters 


this looks more like 


Step 3: evaluate whal |’m used to 


© Another 2-d problem 


Evaluate | F.dr 
C 


a 


where F =(ax+y)i- x | 


and C, is the quarter circle, centre af 
the origin and radius 2, belween 


(2,0) and (0,2) 
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Equation >,* + _" a=: 4 


x = 


ue 
C, starts af (2,0) and 


ends. at (0,2) 
So t goes From 


io f= Eo 


© Substitution 
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(-4-8 sin t cost ) dt 


Solution at the end of the Frames 


® Comparison of answers 


y 
(0,2) Le 


® Another vector field 
evaluate i F.dr 
& 


where F= x71 * ¥) 


a) for C=C, (li) For a 


26 
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met. ye ent 
pe i +o) 
oe. 
si 
C2 a. 
{ E-dr = | (E-dr)at= | ¢ +t-2)dt=-4 
C, C, dt > —_ 
Fa 
diy C=C, 
x=2Zcosl y= est 
E =A costtit 2sin tj 
oe 4 foe | 
at =-2sinti +2 cos t | 
bt /2 
| Fede ome (F )dt = |(-Beoss ct sin +4 cos E sin t) dt = -4 
a : , 3 
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4 A closed path 


@ Fedr= 
Cc 


where F = yt + 2) 


and C is the closed curve 


APBQA 
Two 
noe Bede =[ fide + { F. dr 
paths c APB BQA 
a cs a Sse naam ern 
at) { F. dr : (ii) { Fe 
APB BQA 
Parametrize | Parametrize 
(0,2) | (0.2) 
| 
(2,0) | (2,0) 
_ gubstitute | substitute 
t 
[le #2) de , {*(e. 2) at 
b dt | ; c 
1 ] 
evaluate | evalvafe 
| 
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Answers 


(i) | Fae 
APB 


Porametrize 


y te x «2 cost 
4 y = 2 sin 
a 7 AY 
c — O and Ct, = /> 
x 
substitute 
fF = 4 om CU cost it 2) 
ae aes et : 
qe iN EL +2 cost) 
evalvale 
N/2 4 
{ (-8 sin? Ecos t +4 cos t) See 3 
e) 
(11) | F.dr 
BQA 
Parametrize 
Y Ts x= 
Q y = 2c 
— C = O and. fae 
substitute 
F = (2-0 se 
de 
ac = ae 
evaluate 
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for every 


CONSERVATIVE FIELD 
closed path , 


closed path ) NON-CONSERVATIVE FIELD 


17) Summary : How To evaluate line integrals 


Step | > parametrize 


step 2 : substitute 


slep 5 : evaluate 


This applies to three—dimensional line integrals 


as well as two-dimensional line integrals 
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Solutions to Exercises in the tape 


Frame 5 
F= (24+ ti-tj 
r=ti+ (2—-—1)j 


d 

Pw 24+ Ht] 0+ Ww 
dt 
L_-—t ta a ao 


0 
{ (2 + 2t)dt 
2 


Frame 11 
n/2 n/2 
(—4 — 8sintcost)dt = [41 — sin? 


0 0 


= (—2n — 4) - (0) 
= —4—2rn. 


Frame 15 


The path APB is the same as the path C, and Frame 13a gives the solution to the line 
integral of F = x?i+ yj along C,. We have 


| F-dr = F-dr = —%. 
C2 APB 


The path BQA is the same as the path C, but here we go along the path from B to A. In 


Frame 13a we integrated from A to B and found that | F-dr = F-dr = —#. 
, C1 AQB 


The integral from B to A is just minus this answer because we just reverse the limits. 


Thus F-dr = 4. 
BOA 


3.4 Post-tape exercises 
Exercise 2 


Evaluate the scalar line integral u- dr for each of the vector field functions u given below 
and where C is a circle defined by the parametric equations | 
x =acost, y = asint, z= (0 <= t<. 20). 
fi) w= (ei + 9 
(ii) u = —yi+ Xj. 
[Solution on p. 45 | 


The solution to this exercise illustrates one important point that was introduced 
on the tape: that the value of a line integral round a closed curve may not be zero. 


This exercise and the problems in the tape section are examples of line integrals 
along paths drawn in the x,y-plane — these we called two-dimensional line 
integrals. We shall now look at some paths in three dimensions. The same three- 
step approach that we introduced on the tape applies here; the steps are 

(i) parametrize, (ii) substitute and (111) evaluate. 


Example 3 


Evaluate the scalar line integral of F = xyi + yz*j + y*zk along the path C 
between the points (0,0,0) and (1,1,1) defined by the parametrization 


x = t, y=,. z=t << ts 1). 


MST204 26.3 31 


Solution 


The parametrization of the path is given in the question; this is usually the case 
for curves in three dimensions —- they are defined in terms of a parameter. The 
limits of the integration are t, = 0 and t, = 1. 


Now we must replace x, y and z throughout in terms of t. We have 
F=fi+r+j+ek 

and r=xi+yj+zk=ti+ t?j+tk 

so that 


ie , 
att tit k 


d 
and Fo = 4254 8 = 2 4 345. 


On substitution the line integral | F- dr becomes 
Cc 


1 
| (t? + 32t5)dt. 
0 
This is a familiar looking integral and evaluating it in the usual way we get 


1 
[ Fear= [ (t? + 3t°)dt = 3. 
¢c 0 


Exercise 3 

Evaluate the scalar line integral of 
F = (2x)i + (xz — 2)j + xyk 

along the path C, between the points (0,0,0) and (1, 1, 1) defined by the parametrization 
x=( y=, 2=f (0<t< 1). 


[Solution on p. 45] 


Exercise 4 


In Exercise 3 the path of the integration starts at (0,0,0) and ends at (1,1, 1). Repeat the 
exercise for the path starting at (1,1, 1) and ending at (0,0,0) and following the path of 
Exercise 3 in the opposite direction. 


[Solution on p. 45] 
The solution to Exercise 4 illustrates an important property of the scalar line A 


integral. If we reverse the end points of the path then the value of the line integral 
changes sign, i.e. for the path C in Figure 6 


F-dr= — F: dr. 


AB BA - 


Figure 6 
Exercise 5 
(i) For the vector field defined in Exercise 3, evaluate | F- dr for the path C, between 


C 
the points (0,0,0) and (1,1,1) defined by the parametric equations 
X= yet <gef <t<1). 
(i!) Use the results of Exercise 3 and part (i) to find the value of ‘ F-dr where C is the 


Cc 
closed curve starting at (0,0,0) going along C, to (1,1,1) and returning to (0,0,0) 
along the path C,. 


[Solution on p. 45] 


Example 4 


Show that the scalar line integral F’- dr for the vector field F = x?i + yj is 


AB 
independent of the path taken between A and B. 
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Solution 


Consider two points A and B with coordinates (x,,y,,Z;) and (x2, y2,Z2) 
respectively, and any curve C joining A to B (see Figure 7). 


Let a parametrization of the curve C be 

x=f(t), y=st), , z= h(t), 
with t = t, and t = t, at A and B respectively. 
Then F = (f(t))?i + g(t)j, 


ee dg. dh 


dt Fe at? re, 


itary { (t2))? (f(t1))? . (g(ts))? 
= a + g 5 | _ (eur + = 


Now f(t,) and f(t) are the x-coordinates of the points A and B; also g(t,) and 
g(t) are the y-coordinates of the points A and B. The value of the line integral can 
thus be written in terms of the coordinates of A and B as follows: 


3 2 

y2 Xi Vi 
F- ee) (gears ean’ 
- dr = 343) 43 


This result depends only on the positions of the end points A and B and not on 
the path joining A and B. 


The vector field in Exercise 5 is an example of a non-conservative vector field 
because in part (ii) we found that the value of its line integral round a closed curve 
is non-zero. Compare this with the vector field defined in Example 4. Consider any 
two paths labelled APB and AQB joining A to B (see Figure 8). The value of the 
line integral along each of these two paths is the same: as we found in Example 4 
it is 


Hence the value of the line integral of x7i + yj round the closed curve APBQA is 
zero because 


{ iat = F-dr+ F-dr = F-dr — F-dr = 0, 
APBQA APB BQA APB 


AQB 


and this will be true for every closed path. The vector field F = x?i + yj is an 
example of a conservative vector field. 


Exercise % 


(i) Show that the vector field grad ¢ where ¢ = a fal 


ae = is the same as the vector field 
given in Example S, 
(ii) Show that for the vector field grad ¢ in part (i) and the points A = (x,,y,,z,) and 


= (X2,¥2,Z2), 
| grad b: dr = $(x2,y2,Z2) — 6(%1, 91,21). 
AB 


[Solution on p. 46] 


xX 


Figure 7 


B(x ,¥,Z>) 


x 


Figure 8 
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Exercise 7 


Use the same method as in Example 4 to show that the value of the line integral of the 
vector field grad @ for any scalar field is given by ¢(B) — (A), where A and B are the two 
end points of any path in the domain of @. 


[Solution on p. 46] 


The result of Exercise 7 shows that the vector field grad ¢ obtained from a scalar 
field @ is a conservative vector field. We shall continue the discussion of 
conservative vector fields in Subsection 4.3 of this unit. 


Summary of Section 3 


1. The scalar line integral is of the form | F’- dr where F is a vector field function 
and C is a curve along which the diese is carried out. If C is a closed 
curve then we write F-dr. 

“ C 

2. A curve can be expressed in terms of a parameter t by the parametric 

equations 
x = x(t), y = y(t), z = 2(t) (a<t<b).. 

3. To evaluate a scalar line integral we express F and r in terms of a parameter t 

and then the line integral F- dr can be written in terms of an ordinary 


integral in t, as summarized in the following procedure: 


a eae RS OI 


Procedure 3.4: Evaluating line integrals 


Step 1. Parametrize: find a parametrization of the curve so that 
every point on the curve is given in terms of a parameter 
3 


Step 2. Substitute: write each variable x, y, z, the vector field F 
dr 
dt 
an integral in terms of t. 


and — in terms of t and write the line integral | F + dr as 


Step 3. Evaluate: evaluate the t integral. 
ee 


4. The line integral of the gradient of a scalar field ¢ is independent of the path. 
We have Vo: dr = $(B) — (A). 
B 


A 


S. A vector field F is conservative if for every closed path C drawn in the domain 
of F we have 


} F-dr=0 
| & 


Otherwise F is a non-conservative field. 


4 The curl of a vector field 


4.1 Line integrals and curl (Television Subsection) 


The television programme investigates line integrals of velocity fields representing 
three fluid flows and uses the line integral of a vector field u to define a new vector 
field called curl u. 


During the television programme the results of the following exercises are used; so 
work through the exercises before watching the programme. 
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Exercise 1 

A circular disc rotates about its centre O with constant angular velocity w; that is to say, 
the angle turned through in time t is wt (radians). Any point P on the disc has velocity of 
magnitude wa | and in a direction perpendicular to the radius vector OP where a is the 
magnitude of OP (see Figure 1). 


(i) Show that the velocity vector field u is given by u = w(—yi-+ xj) using Cartesian 
coordinates with origin at O. 


(ii) Evaluate the line integral u-dr where C is a circle whose centre is the origin and O x 
Cc . 
whose radius is b. Figure 1 


[Solution on p. 46 | 


Exercise 2 : 
If f(y) is a differentiable function of y, show, by using Taylor series, that 


i dy\ | Pee dy If you find this exercise 
i ag 2 df difficult, consult the solution. 
ee ee ee 

dy70 oy dy 


[Solution on p. 46 | 


The programme can be summarized as follows. 
Part 1: we begin by introducing two different flows of a liquid: 


(i) the motion of water in a bowl which is rotating with constant angular 
velocity, and 


(ii) the motion of water in a vortex tank, which is rather like what happens to 
bath water near to a plughole. 


These two flows appear similar because the motion is circular; but small floats 
placed in the fluid behave differently in the two cases. 


Part 2: the rotation of a rigid disc is used as a simple model for the motion in the 
bowl and we find the velocity field u representing this flow. We then use a line 
integral to define the z component of a new vector field called curlu. The 
expression defining (curlu), relates the line integral to the rotational properties of 
the vector field u. The value of (curlu), is used to predict whether small floats 
placed in the fluid flows will rotate or not. 


Part 3: the flow in a long rectangular tank is then investigated. At first sight there 
appears to be no rotational motion of the fluid in the tank and we would expect a 
float not to rotate. However the formula for (curlu), predicts otherwise. Have we 
made a mistake? 


Now watch the television programme and find out. 


In tke television programme we defined the z component of the curl of a two- 
dimensional vector field in the x,y-plane. We can define the curl for any such 
vector field in the following way: | 


TV26 


Suppose that F is a two-dimensional vector field defined over a 
region of the x,y-plane. Then the z component of the vector field 
curlF at a point P in the domain of F, is defined by 


1 
(curlF), = lim (<4 Far 
A~0\AJdc 


where C is a curve drawn in the domain of F surrounding the 


point P and A is the area within C (see Figure 2). 
. | 
The result is the same regardless of the shape of the curve C. 3 


an On nan 


Figure 2 


During the programme the results of several calculations involving the limit in this 
definition above were assumed. In the following exercises you are asked to show 
that these assumptions are correct. | 
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Exercise 3 
The velocity field for a disc rotating about the origin with constant angular velocity w is 


1 
given by u= — wyi + wxj. Show that lim (<4 n° | = 2w in each of the following cases: 
A-0O C 


(i) C is a square of side b, and centred at the origin (see Figure 3); 
(ii) C is a square of side b, and centred at a point (xo, yo) (see Figure 4). 


[Solution on p. 46] 


Exercise 4 
The velocity vector field for the flow in the vortex tank is given by 
ts 


rae 


where e@, is a unit vector perpendicular to the radius vector at a point P (see Figure 5) . 
distance d from the origin. 


—yi+t+ xj 
Jx? + y? 


(11) Hence show that u = 


(i) Show that e, = 
c 

eo (-—yi +x j). 

(iii) Evaluate the line integral of u round the curve C = ABCD shown in Figure 6, where 
AB and CD are parts of radii and BC and DA are arcs of circles. 

(iv) Hence evaluate (curlu), at the point P. 


[Solution on p. 47] 


We saw in the television programme that the curl of a velocity field u describes the 
rotational properties of the field. We could predict whether a small float placed in 
a fluid would rotate or not depending on whether the value of (curlu), was non- 
zero or zero respectively. 


The mathematical analysis for the third fluid experiment, involving the linear tank, 
introduced a different way of calculating curl. We found that for the velocity field 


d 
u = f(y)i, the z component of curlu could be expressed as ->. 1.e. in terms of the 


rate of change of the velocity field. 


Calculating derivatives is certainly easier than evaluating line integrals and taking 
limits. In the following example we obtain a formula for (curlF), for the more 
general two-dimensional vector field F = F,(x, y)i + F,(x, y)j, in terms of the 
spatial rates of change of the components of F. 


Example 1 


Suppose that F is a two-dimensional vector field defined over a region of the x, y- 
plane. By integrating round a small rectangle show that . 


A 
a 
Figure 3 


D & 


> 


O x 
Figure 4 


O 
Figure 5 


O 
Figure 6 
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“4 
oy 


where F, and F, are the x and y components of the vector field. 


ao see 


Solution 


To find the z component of curlF we need to evaluate lim 36 F- dr for some 
Aw0A 

area A contained within a closed curve C. Consider a two-dimensional vector field 
F defined over a region of the x,y-plane and let F = F,i + F,j. For the curve C, 
consider a small rectangle ABCD with sides h and | and centre at the point P with 
coordinates (Xo, Yo) (see Figure 7). To evaluate the line integral of F round the 
rectangle we shall appeal to the definition of a line integral as the sum of terms of 
the form F: or where or is a small displacement along the curve of integration. For Figure 7 
the rectangle ABCD, since we are assuming that h and Tare small, suppose that we 


divide the rectangle into four small displacements AB, BC, CD, DA where 


AB = hi, BC = 1j, CD = —hi and DA = —lj. 


We now approximate the line integral as the sum of four scalar products in the 
following way: 


F-dr ~F,,°AB + KF, BC 4 F.C) + F.,-DA 
Cc 


where F AB> F gc, F cp and Fp, are the values of F at the mid-points of the sides of 
the rectangle. | 


Now 
ey l 
F 453° AB =F ,,°hi = Filxo.y = )h 


because F ,,°i is the x-component of F evaluated at the mid-point of AB, which 


l 
has coordinates x, Yo ). 
aoe a h 
Similarly F gc-°BC = F(x + 5-99} 
pokes ] 
Fep° CD — ~Fi(xo.¥0 + 5}) ee 


pos h 
Fp,* DA = ~Fi{xe a + yo} | 


and so we have | 


| 
oF dr = [Fi(x0.90 - )- Filxo.r0+ 5})h 
h 
+ [Blo 590) — Flo - > vo} 


The area contained within C is the area of the rectangle ABCD and this is hl. Thus 


, l l 
Fi(xo.r = 4 me Filo + 5 


1 f 
pF dr = ; 


h h 
Fx. = i 5-99 = A Xo — 599 
7 h 


Now in the limit as A ——> 0, both h—- 0 and !——>. 0; and we have 


aa 
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l 
Fi(xo.vo — 3] ~F.{ sore + 3 OF 
= 1 


h,l>0 oy These results follow either by 
expanding the numerator as a 
Taylor series about (xo, yo) 
and and then taking the limit as 
h h Y Z <0), v2 has 
ee see ese 3 ennition of partia 
F[x 3 ad | F[x “J ; OF, differentiation in Unit 25. 


lim 
h,l-O 


h 
1 OF, OF 
Finally we can write lim apr ‘dr = ——~ + —* and so 


A~>0 


OY 5 Ox 


for a two-dimensional vector field. 


OF 
Note that if F, = 0 then (curlF), = aay which is in agreement with the result 


obtained in the television programme for the velocity field of the flow in the linear 
tank. 


We now have a formula for the z component of curlF in terms of the derivatives of 
the components of F. In the next exercise you can check the results of Exercises 3 
and 4 by using this derivative formula. 


Exercise 5 
Use the derivative formula to find the value of (curlu), for the velocity vector fields 


(i) u= —oyi+ oxj; 


= c . , 

i) u=—.——_(_— ypi+ xj). 

(11) eet yy! y j) 

In each case compare your answer with the relevant result obtained in Exercise 3 or 4. 


[Solution on p. 47] 


The velocity field given in the first part of this exercise represents the motion of 
the rotating disc. We find that (curlu), = 2w at every point on the disc, and that is 
the result we obtained in the television programme by using line integrals. The 
velocity field given in the second part of this exercise represents the flow in the 
vortex tank. We find that (curl u), = 0 in this case, which again agrees with the 
line integral calculation. 


So the curl of a vector field describes the rotational properties of the vector field. 
Its effect can be illustrated by placing small floats into a moving liquid. Then the 
curl of the velocity field at a particular point is twice the angular velocity of a float 
placed at that point. 


4.2 Curl of a three-dimensional vector field 


The television programme and the definition of curl in the last subsection 
considered only two-dimensional vector fields defined over regions of the x,y-plane. 
We defined the z component of curl for a vector field F = F,i+ Fj by the line 
integral 


lim —Q F:dr 
AO C x 
Now the z direction is normal to the area A (see Figure 8); and the unit vector k Figure 8 


points in the z direction so we could have written (curlF), as (curlF)-k. 


We define curl for a three-dimensional vector field in an analogous way. 
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Suppose that F is a three-dimensional vector field defined over a 
region of space. Then there exists a vector field curlF with the 
following property: for any point P and any unit vector n, the 
component of curlF in the direction of n is given by 


(curlF), = (curlF)-n = lim@® F-dr 
A707 C 
where C is a small curve lying in a plane and enclosing the point 
P such that n is normal to this plane containing C and A is the 
area of the plane enclosed within C (see Figure 9). 


This is rather a formal definition which is a simple extension of the definition of 
curl for a two-dimensional vector field. You are more familiar with vector fields 
being given in terms of a vector function and its domain. The domain of curlF is 
the same as the domain of F because curl is defined for the same set of points as is 
F. To write curlF as a vector function of the form a,i+ a,j + a;k we need the 
components of curl in the directions of i, j and k. Now we can find these 
components by using the definition above: 


1 
a, = (curlF)-i= lim —@ F-dr, 
A10A1 Jc, 


1 
a, = (curlF)-j = lim —@ Fv-dr 


A270A2 Jc, 
1 
and a; = (curlF)-k = lim —@ Fv-dr 
A370 443703 


where C,, C, and C; are curves drawn parallel to the y,z-plane, x,z-plane and x,y- 
plane respectively, and A,, A, and A; are the areas contained within C,, Cz and 
C,; respectively (see Figure 10). 


Figure 10 z 


x 


So we can write 


ee | ee. 
curl = { im 7 Far + {tim 7 Far) 
Ai7O0 444 7 C, A270 442 7 C2 
og 
+ {im F dr) k 
A370 443 ¢ C3 
For a two-dimensional vector field we found that we could write the z component 
of curl in terms of derivatives. Now we show that we can do the same here. 


Example 2 
By taking a rectangle parallel to the x,y-plane show that 
1 OF, OF, 
lim —® F-dr =— —-—— 
A377 r@) A; C3 Ox oy 


where F, and F; are the i and j components of F respectively. 


Solution 
Consider a vector field F and let 


F =F,i+ F,j+ Bk. 


x 


Figure 9 


When n is the unit vector k, 
then (curlF), is just what we 
denoted previously by 

(curl F),. 
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For the curve C3 take a small rectangle ABCD with sides h and / and centre at the 
point P with coordinates (xo, yo,Zo), the rectangle being drawn in the pees 
Z = Zo = constant (see Figure 11). 


The analysis is similar to that used in the solution to Example 1. We write the line 


integral F-dr as the sum of four scalar products as follows: 


x 


$ F-dr ~F 43: AB + F pc: BC + Fcp*CD + Fp,*DA. Figure 11 
C3 


For this vector field F we have three components; however the four scalar 
products involve multiplying by either i or j so that the k component of F, ie. F;, 
does not appear in the value of the line integral. For example, 


ices ] 
F43° AB = Fy 3° (hi) = Filxoon 5 520)h 


Sa h 
and Fgc* BC =F 3° (lj) = Fl + 5 Yonzo]l 


| | 
1 Filxoon ~ 5.20] ~ F.{ sanyo + 3420 
Tht — t= ..L.". e 


3¢7 C3 l 
h h 
Fl = 5 Ponto 3 Flo “ 5 Yost 
oa 
h 
and in the limit as A; 0 (i.e. h——0 and 1!——>.0) we have 
1 OF, OF, 

lim —® F-dr =—* —- — 

4370 A3 cx. Oy 
and this is the k component of curlF. In a similar way we can obtain derivative x 
forms for the i and j components of curlF. For the i component we choose for C;, Figure 12 
a small rectangle in the x = xq plane (see Figure 12). The i component becomes 
OF, OF. , . 
rr ~ oo For the } component we choose for C, a small rectangle in the y = yo 

y Z 
OF, OF 
plane (see Figure 13). The j component becomes a — = 
Thus we can write the vector field curl F in a form that is easy to evaluate, namely: 
x 
Figure 13 


OF, OF, = OF, OF; ~ OF, OF, 
curl F = —]i 


oy = az éz oxi? \ax dy 


Applying this formula to the two-dimensional vector field F = F,(x, y)i + F,(x, y)j 


OF, OF 
we obtain the same result as in the last subsection; that is (curlF), = ae ae 
but now we can go one stage further and write curlF as the vector 

0 0 
curlF = = Se ll 
oy 


We were not able to do this before because we knew nothing about the x and y 
components of curlF. However we could have made an educated guess at their 
values. We saw in the television programme that for the rotating disc, 

(curlu), = 2m. Now angular velocity is a vector quantity; for the rotation of the 
x, y-plane the angular velocity vector is wk, so we might have suspected that 
curlu is also in the k direction. 
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Example 3 

If F = (xy + 27)i + x?j + (xz — 2)k 
(i) find curlF at the point (1, —1, 3); 
(ii) find points where curlF = 0. 


Solution 


We use the differential formula to calculate curlF. 


F,=xy+2, £¥F, =x’, F, = xz — 2, 


0 
so that of, = x and * es 
oy of 
oF. OF, 
3, ~ x and = 9, | 
OF, OF; 
Ae = 27 and Gy: = 0, 


giving curlF = (0)i + (2z — z)j + (2x — x)k. 
(i) At the point (1, —1,3) we have curlF = 3j +k. 


(1) curlF is zero when x = 0 and z = 0 for any value of y. 


Exercise 6 

Find curlF if 

(i) F=r; 

Gi) F =Jirjt; 

(iii) F = (y? + 2z)i + (xy + 6z)j + (2? + 2xz + y)k; 

(iv) F=grad@_ where ¢ is a scalar field. 

Herer=xi+ yj+zk,r=[r|= Jets and f is a unit vector in the direction of r. 
[Solution on p. 48 | 


The type of vector field function in part (ii) of this exercise can be used to 
represent many physical quantities in nature. For example, you will see in a later 


GM 
unit that the gravitational force field can be written in the form F, = — r ft SO 
GM : ; 
that f(r) = — steed The result of this exercise shows that curlF, = 0. 
r 


We will see in the next subsection that if the curl of a vector field is zero for all 
points in its domain, then the vector field is conservative, so we will use curl 
mainly as a test for conservative fields. 


The derivative formula for curlF contains terms involving partial derivatives. The 
expression is not particularly easy to remember; however the following exercise 
shows an alternative notation in terms of the differential operator del. 


Exercise 7 
By expanding the cross product of V and F show that curlF = V x F. 
[Solution on p. 48 | 


This is an alternative definition for the curl of a vector field which you may come 
across in other texts. The way we have defined the curl in this section more readily 
shows its important physical significance. 


4.3 Conservative fields 


C(xi+ yj +k) . . 
Se is one that we have considered at various 
hae ae eee 


times in this unit. In Example 2 of Section 1 we saw that this vector field function 


The vector field F = 
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represents the gravitational force on a body a distance r from the Earth’s centre. 
In Example 3 of Section 2 we showed that it could be derived from the scalar field 


C ; = 
o = ————-. In Exercise 6 of Section 3 we saw that the line integral of F 
./ x? + y? + 2? 


between two points A and B is independent of the path taken. And Exercise 5(ii) 
of this section shows that curlF = 0. Now we shall put these ideas together by 
talking about a special class of vector fields called conservative fields for which 
these last three properties are true. 


On the audio-tape we defined a conservative vector field F as one for which the line 
integral F'- dr = 0 for every closed curve C drawn in the domain of F. (The 
c 


reason for the name has to do with the conservation of energy for bodies acted on 
by forces which are modelled by conservative vector fields.) If a vector field does 
not have this property then it is non-conservative. The velocity field u representing 
the rotating disc is an example of a non-conservative vector field, because when C 
is a circle of radius b and centre the origin, 


4 u-dr = 2wnb? (see Exercise 1). 
Jc 


This definition of a conservative vector field has an important consequence for the 
curl of F. We have seen that 


1 
curlF = lim [7-4 Far) + lim (+4 Far) 
Ai;7>0 A, Cy A2—>0 A). cs 


1 
+ lim (7-4 Far k. 
A3>0 A; C3 
Now for a conservative vector field, the line integral round every closed curve is 


zero. This means that each of the components of curlF is zero, and we have the 
important result: 


If F is a conservative vector field then curlF = 0. 


The converse of this is also true: if curlF = 0 and the domain of F is sufficiently We shall not say precisely 

simple then F is conservative; however we will not be able to prove this statement what we mean by ‘sufficiently 

aoe ee simple’. Certainly every 

in this course. This gives an alternative definition of a conservative vector field domain we conquer i dik 

which you may come across in other texts. unit is sufficiently simple in 
this context. 


lf curlF = 0 at every point of the domain of F, then F is said to 
be a conservative vector field. 


Thus if we want to decide whether a vector field F is conservative or not, instead 
of evaluating line integrals round every possible closed curve that could be drawn 
in the domain of the vector field, we evaluate curlF in its derivative form and see if 
its value is zero or not. 


Example 4 
Which of the following vector fields is conservative? 


(i) F = zx*i+ zy?j +4(x? + y’)k; 
(ii) G=xit x’yj. 
Solution | 
We have to evaluate curlF and curlG. Using the derivative formula for curl gives: 
(i) curlF = (y* — y*)i + (x* — x?)j + (0 — 0)k 
= 01+ 0j + 0k = 0. 


Thus curlF = 0 so that F is a conservative vector field: 
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(i) curlG = (0 — 0)i + (0 — 0)j + (2xy — 1)k 
= (2xy — 1)k. 


Thus curlG # 0 at every point, so G is a non-conservative vector field. 


Exercise 8 


(i) The vector field F = 2xi + (xz — 2)j + xyk was shown to be non-conservative in 
Section 3 (Exercise 3) by evaluating the line integral along a closed path. Check this 
result by finding the value of curlF. 


(ii) Is the vector field G =r"f conservative, where r = xi+ yj+zk, r= |r| and f =r/r? 


[Solution on p. 48 | 


Exercise 9 
If F is a conservative vector field show that the line integral of F between two points A and 


B, 1; [ F : dr, is independent of the path taken. 


* AB 
[Solution on p. 48 | 
Exercise 9 shows that for a conservative vector field the value of the line 
integral a dr depends only on the end points A and B of the path and not on 


- path itself. In particular we can define a scalar field ¢ at every point 
= (x,y,z) in the domain of the conservative field F by letting 


wx 2) = | F< dr (see Figure 14). 
OP 


XxX 


(1) Figure 14 


Then | — = $(B) — @(A). 


We recall from the result of Exercise 7 in Section 3 that the value of the line 
integral of grad @ along any curve between two points A and B is @(B) — @(A). 
That is, 


{ grad d-dr = ¢(B) — (A). 

AB 

In general it can be shown that when a (1) hold then F = grad @ for the 
scalar field @. 


One application of this idea is to the motion of a particle acted upon by a 
position-dependent force which can be modelled by a conservative vector field. 
Such a vector field is called a conservative force field. For instance, consider a 
body moving in a one-dimensional force field F = F(x)i. Any one-dimensional 
force field is conservative and the scalar field which we define from this force field 


is just d(x) = | F(s)ds. In Unit 4 we defined the potential energy of a body to be 
0 


U(x)= —- : F(s)ds, so the one-dimensional scalar field that we have introduced 
here is just ee the potential energy of mechanics. 

We can generalize the ideas of one-dimensional mechanics and we use the line 
integral — [F -dr to be the potential energy U of a body moving in a conservative 


force field F. Conversely, given a potential energy function U(x, y,z), the force 
acting on a particle is 


F = —VU. 


In a later unit we will follow this idea up. You will see that for a conservative force 
field we have a conservation of energy equation, 


kinetic energy + potential energy = constant, 


and thus the name ‘conservative’. 
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Summary of Section 4 


2 


For a two-dimensional vector field F defined over a region of the x,y-plane, 
the z component of curlF at a point P is defined by 


1 
(curlF). = lim Fr Far 
A-0 A eS 


where C is a small curve drawn in the domain of F surrounding the point P 
and A is the area within C. If F = Fji+ F,j then an alternative way of 
writing (curlF), is 


For a three-dimensional vector field F we define the component of curlF at a 
point P in the direction of a unit vector n to be 


17 : 
(curlF), = (curlF)-n = lim “4 Far 
A-0 A ‘a 


where C is a small planar curve enclosing the point P, A is the area enclosed 
within C and n is normal to the plane of the curve C. 


We evaluate curl of a vector field by using the derivative formula: if 
F=F,i+ F,j+ Bk 


then curlF = V x F = A). OF, OF). ee A 
fz. Gx Ox. oy 


oy Oz 
The curl of a vector field F is another vector field defined over the same 
domain as F; curl is a detector for rotational motion. F or instance, the curl of 
a fluid velocity field describes the spin of a float placed in the fluid. 


For a conservative vector field: 


(i) curlF=0 (and conversely, if the domain of F is sufficiently simple); 
(11) | F- dr = $(B) — $(A) where F = grad ¢. 
AB 


The curl of a vector field F can be used as a test to see if F is conservative. 


44 MST204 26 Solutions 


Appendix: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. The temperature T at a point r = xi+ yj on the metal velocity vector field u is directed along OP and has 
sheet is defined by fges log. r. But T = 0 when r = 1, and magnitude c/loP| for some constant C. So 
T = 100 when r = 4. Putting T = a(log,r) — b we get b = 0 pies 
and a = 100/log.4. Thus the scalar field function is COP C 
u(x, y) = = = Pas eee Te i+ yj). 
100 log,r loP||lop| +’) 
fags. 5. The Earth’s gravitational force field is directed towards 


See the centre of the earth. The left-hand figure shows arrows” 
Since r = ,/x* + y”, the scalar field is defined over the region _ representing the gravitational force at various points. The 
of the plane given by the set of points (x, y) such that field lines are shown in the right-hand figure. The 


l<./x + = 4. 7 arrowheads indicate the direction of the force field. 
2. The contour curves are given by 


xy = constant =c 


c 
le. y =—. 
x 


6. If you look back at Exercise 4 the field lines in this 
exercise would describe the velocity vector field of that 
exercise. So the direction of the vector at each point is in the 
direction of xi + yj. The vector field function has the form 
F = f(x, y)(xi + yj) where f(x, y) is some function of x and y 
taking positive values only. 


Solutions to the exercises in Section 2 


The figure shows the contour curves for c = 4, c = —% and 
C= t. 1. ) << ys 
3. (i)(a) There is a hill on the left, near ‘contour 100’. ad wr) 

—=2xy and = —4 at (-1,2); 
(b) There may be a lake where the contours are marked with Ox ax 
negative numbers indicating land levels beneath sea level. ad ad 

aa eae ee 

(ii)The land should be fairly level on the right; the contour oy sie and oy 1 at (—1,2). 
curves join places at levels 0, 1 and 2 over quite a large area. So 
(iii) If we start off at point A and walk in the direction of the 
arrow the land begins to climb very quickly; the contours are grad = se 4 aa = —dj + j at (—1,2} 
close together so the heights change very quickly indicating a oy 
hard climb! 


2. The curve x” — 2xy + y* = 9 is a contour of the scalar 
4. The vector quantity is the velocity of the water at points field @ = x” — 2xy + y?. The vector grad ¢ evaluated at (0, 3) 


on the plate. The vector field function u (i.e. the velocity is normal to this scalar field; 
field) is defined for points P on the plate. The region over = : ; 
which the vector field is defined is the set of points (x, y) eee ex es 
inside the circle x? + y” = 100 and excluding the origin. The and at (0,3), grad@ = —6i + 6j, and this is the vector 
required. 
Ei 
a2 b= xy + 2X; 
0p 0d 
| 2 he ee = 
ax xy +z and ax 2 at (—1,1,0); 
ao _ hos 
ee d Seas oe ; 
By * an By 1 at (—1,1,0); 
0g oy) 


— =2xz and —=0 at (—1, 1,0). 


MST204 26 Solutions 


So grad @ = —2i1+j + Ok at (—1, 1,0). 
4. (i) grad db = 2xy?z?i + 2x*yz? j + 2x?y2zk. 


The maximum rate of change is |grad @| in the direction of 
grad d 


lgrad g| 
lgrad d| = |—2i1 + 2) + 2k| =2,/3 


At the point (—1,1,1) we have 


and 
grad @ Se 1 
ears = Wes i We 
(ii) The rate of change of ¢ in the direction of e is (grad ¢)-e. 


For e = rt a 5* we have 


6 8 
(grad d)-e = gxyez - gx "2. 


At the point (2,1, —1) we have (grad d)-e = —4. 
5. For two scalar fields ¢ and w we would expect that 


V(oy) = (Vow + O(V). 


To check this result we consider V(dw) in terms of the 
partial derivatives. 


6 0 7) 
V(ow) = = (oW)i + — (GW)j + — (OW)k 
Ox oy Oz 


_ (6¢ ow\. | [od oy). 
= (av rea) (S405) 


_ 


og oy 
an + 62 k 


Op, ow, ow 
+ fein P+ Me 


= (Voy + O(V). 


Solutions to the exercises in Section 3 


| ee (a,0,0); 
t = 7/4, (a/,/2, a/,/2,0); 
t = 2/2 (0, a, 0); 
t = 3n/4, (—a/,/2, a/./2,0); 
t=, (—a,0,0); 
t = 3n/2, (0, —a,0). 


The curve lies in the plane z = 0, and is a circle of radius a 
with centre at the origin. 
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2 Gi) x=<acost, y =asint; 


tL =, ty = 28: 


1 e . e 
u= q2acos ti + asintj); 
— = —asinti+ acostj; 


f u:dr = oa z= Q. 


“C 0 

(ii) The parametrization is the same, so now 
u = —asinti+ acostj; 
dr 


— = —asinti+acostj; 
dt ee: ; 


2n 
u-dr = | a*dt = 2na?. 
C 0 


3. x=t, per, z=P; 


i, = 6, L_=f. 


F = 2ti + (t* — 2)j + ek; 


d 
= =i + 2yj + 30k. 


dt 
1 
F-dr = [ (St? — 2t)dt = —2. 
Ci 0 
4. x=, y=P, z=: 
t, = 1, t, = 0; 


F = 2ti + (t* — 2)j +k; 


dr 
—=j+ 2tj 303k; 
e? 2+ Zt + 


0 <4 
[F-a= | (5t° — 2t)dt = +2. 
1 


5. C, is defined in Exercise 3 


Zz 


x 


(i) Line integral along C,: 


x=t . y=t. z=’; 


t, = 0, t, =1; 
F = 2ti+ (t°? — 2)j +-27k; 


dr 
—=i1+j+2tk; 
i 26 Bae 


1 
F-dr = } (2t + 3t° — 2)dt = —4. 
C2 0 
(ii) For the curve C, we go along C, and return along C,, so 


F-dr = —¢—(—2) = x. 


C2 


b F-dr = F-dr — 
C Cy 


46 
e x" rs 
° If = — ees 
6 @d=>+5 
7) 
then grad ¢ = aa + @ + ae 
oy Oz 
= x7it+ yj 


and this is the vector field used in Example 4. 


(ii) From Example 4, 
3 3 2 
x 93 “og eS 
F-dr = |— +=] -|-+> 
heals 45] (5%a) 
oe P(X2,2,2Z2) = P(X1,Y1, 21): 
7. Suppose that A and B are any two points in the domain 
of a scalar field @ and let C be any curve joining A to B. 


Zz 


B 
C 
y 
x A 
Let a parametrization of the curve be 
x=f), y=st, z= h(t), 
with t = t,; and t = t, at A and B respectively. 
0p.  o¢. oi 
Then grad @ = xt + ay! - Ao 
dr dx, dy, 2 
—j+—k 
dt dt a Foe es 
and hence 
grad Rsk Ogdx <dddy <dddz 


ho ba Gee bet 


d 
Using the chain rule from Unit 25, this can be written as -. 


thus 


t> d 
grad d- dr = { era oF dt 
AB dt 


wae : a 
= (tz) — (ts) 
= $(B) — (A). 


This result is true for all paths joining A and B. Hence, the 
value of the line integral of grad @ depends only on the end 
points of the path (i.e. A and B) and not on the path itself. 


Solutions to the exercises in Section 4 


1. (i) The velocity vector at the point P is given in terms of 
i and j by u = —qwasin@i + wacos6j 


and acos@ = x 
asin@ = y- 


so that u = —q@yi + wxj. 
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P(x,y) 


(ii) x = bcost, y = bsint; 
{, = 0, {> = 2a: 


u = —wbsinti+ whcostj; 


— = —bsinti-+ bcostj. 


2 
{ (wb? sin? t + wh? cos? t)dt 


$ u-dr = 
c 0 


2n 
= { wb? dt 
0 


= 2nwb’. 
jy) oy _, dy\? 
2 ve} =soreB rors (2) 2 
dy\* f(y) 
(2 ogo" 


5 2 
and s{y- 2) - so -Zr'0y+ (2) 
(2720, 


5 5 
sly - 3] se a 
ee es 


fr-4) sda 
ad he = fe 


dy0 d 


dy” fy) 
a 


Hence 


3. (i) For the square given in Figure 3, 


$ u-dr = { u-dr+ | u-dr + | u-dr + | u-:dr. 
c AB BC CD DA 


b 
since y= —> 


so taking x as a parameter along here gives 


x =b/2 d b/2 b b2 
| u-dr = | {Fr - dx=[ sey See 
AB x=—b/2 dx bj2 2 2 


b 
Along AB, r = xi — 5) 
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: 
parameter along here, giving 


/2 d b/2 b b2 
| u-dr = [ [r-H}ay = f si el 
BC —b/2 dy —bj2 Z 


b? 
Similarly [ i ad | u-dr = 
2 = 2 


CD 


b b 
Along BC, r=~=i+ yj since x = 5 and y will be a 


Hence } ‘dr = 2wb’. Now area A = b? for the square, so 
im ao ‘dr = B2 + (cob?) se Den: 
ao 


Gi) ude = | Sis vies { u-dr + { u-dr. 
C AB BC CD DA 


b\. 
Along AB, r = xi+ (. -3)i sO 


xo + b/2 b 
| u-dr = | of 0 — 5) a 
AB Xo — b/2 2 
b 
= - 08 — 5. 
a 
Along BC, r= ( + )i + yj, so 
yo+b/2 b b 
| u-dr = [ os +5} dr = 00( x0 +2) 
BC yo- b/2 2 2 


: b\. 
Along CD, r = xi+ (. a s)i Ye) 


xo — b/2 b b 
| u-dr = | ~ ef + 5]ax = 06( 0+ 2). 
CD xy + b/2 2 2 


b 
Along DA, r = (x ~ 5) + yj, so 


yo — b/2 b b 
i u:dr = | ofxo— Say = -06(x9 ~ 3} 
DA yot b/2 2 ps 


1 
Thus u-dr = 2wb? so that lim \ te -dr = 2 again. 
C C 


A-0 


4. (i) e, = —sin6i + cos6j 


and cos@ = aaa sin @ = y 


Ir| Ir| 


x 
Hence e, = ig = 


Moo ee ay 
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(ii) u = =e SO using the result in part (i) we obtain 
un {2+ 3] (where d = ./x? + y’*) 


fo: : 


(ili) Suppose that AD and BC are arcs of circles of radii a 
and b respectively. 


Rg 


p u-dr = [ ude + | u-dr + | u-dr+ | u-dr 
ce AB BC CD DA 


d 
Along AB, u and = are perpendicular so that 


ao 6 sO | u-dr = 0. 
dt AB 


Similarly for the side CD, | u-dr = 0. 


CD 


Along BC, we can use the parametrization 
x = bcost, y = bsint 


and t; = ¢,t,=¢d+4. 
On substitution u = =(—bsin ti + bcostj), 


dr 


— = —bsinti + bcos ti, 
7 sin ti COs tj 


dr ota 
so that ae and j (- Fy at = ca 


Similarly, along DA, 
x = acost, y=asint 


and t; = ¢+4,t, = ¢;s0 


Henoe u-dr =0+ca+0-—ca=0. 
& 


1 
(iv) (curlu), = lim x u-dr = 0 since ucdr = 
A Cc 


A~0 
oF, = ¢ 
5. (curls), = — -— ots 
i oP 
(i) u= —@yi+ wxj 
OF, OF, 
sO ee w and ay = —@ 


so (curl u), = 2a. 


As 


(ii) u Seo: + a) SO 
er re = cox 4 ee 
x (ey) Oty eae 
OF, —¢ 2cy" (y? — x7)c 


oy ty) @+y¥? of +y¥P 
so (curlu), = 0. 
6. (i) F=r=xi+ yj+zk; 

F, =x, Fi, =y, Fy, =2Z; 

curlF = 0. 
(ii) F = f(r; 

=f, B=fO2, K=fo- 
'g r F 

i component of curlF is 


OF; OF, _ - fA 


oy oe ae 


ay ar 


(using the Chain Rule) 


_ BU Se 
“ar r/r var r}r 
Similarly the j and k components of curlF are zero, so 


curl (f(r)f) = 0 for all functions f of r. 


(iii) F = (y? + 2z)i + (xy + 6z)j + (2? + 2xz + y)k; 


OF, OF 
i component of curlF is eae =1-—6= —5: 
oy Oz 
OF, OF 
j component of curlF is — — —* =2—27 = 21 — 2); 
Oz Ox 
7) OF; 
k component of curlF is — y-—2y=-—y; 
Ox oy 
so curlF = —5i+ 2(1 — z)j — yk. 
0 0 
(iv) F = gradgd = - + Be + oe .. 
Oy Oz 
OF; OF, 0 [0p 0 (0p\ | 0: 
dy az ps éz] dz\ay}— 
OF, OF 0 [op op 


@z dx dz\ax -5 es 

OF, OF _ 2 (a¢\ 4 (89) _,, 

dx oy -= éy = 
hence curl (grad @) = 0. 


0 0 0 
7 VxF={i—+j-—+k—] x (Ki+ Fj + Bk) 
Ox oy Oz 
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and using the formula for the cross product given in Unit 14 
this becomes 


(OF: ~ (ah 06) | | (ah 


Siti ee 


8. (i) F = 2xi+ (xz —2)j) + xyk; 
curlF = i(x — x) + j(0 — y) + k(z — 0) 
= —yj + zk. 
Since curlF # 0 in general, F is non-conservative. 


(ii) From Exercise 6(ii) curl(f(r)F) = 0 for all functions f(r). 
For G = rf, f(r) =r" so that curl G = 0 and hence G is 
conservative 


9. If F is conservative then ¢ F-dr = 0 for all closed curves 


zs 
C. Consider a closed curve C passing through the two points 
A and B. 


$ F-ar= | F-de + | F-dr = 0. 
JC APB BQA 
Thus F-dr = ~ { F - dr. 

APB BQA 


But F-dr = — F -dr; 


BQA AQB 


hence F-dr = F -dr and this is true for any path 
APB AQB 


between A and B. So the line integral | F -dr is 


AB: 
independent of the path, and depends only on the end points 
A and B. 


linet Oe TT) 


i 
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